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1 INTRODUCTION 



The first year of observations of the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite has provided us with a 
data set of unprecedented accuracy (Bennett et al. 2003). 
It is exciting to see that we are still far from understand- 
ing all the information that is encoded in the data. Among 
the many papers that have been published based on anal- 
yses of the WMAP data, a large section of them points 
to possible deviations from the standard model. Different 
so called anomalies have been detected and possible ex- 
planations have been offered by several authors: (1) Low 
value of the quadrupole, its alignment with the octopole and 
other multipole alignments (De Oliveira-Costa et al. 2004, 
Schwarz et al. 2004, Vale 2005, Land & Magueijo 2005a), 
(2) Asymmetries (Park 2004, Eriksen et al. 2004a, Eriksen 
et al. 2004b, Larson & Wandelt 2004, Hansen et al. 2004a, 
Hansen et al. 2004b, Donoghue & Donoghue 2005, Tojeiro et 



ABSTRACT 

Motivated by the success of the Bianchi Vllh model in addressing many of the anoma- 
lies observed in the WMAP data (Jaffe et al.), we present calculations in real and in 
wavelet space of the Higher Criticism statistic of the Bianchi corrected Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) first year data. At the wavelet scale of 5 degrees 
the Higher Criticism of the WMAP map drops from a value above the 99% c.l. to a 
value below the 68% c.l. when corrected by the Bianchi template. An important prop- 
erty of the Higher Criticism statistic is its ability to locate the pixels that account for 
the deviation from Gaussianity. The analysis of the uncorrected WMAP data pointed 
to a cold spot in the southern hemisphere centered at I ~ 209°, b ~ —57°. The Higher 
Criticism of the Bianchi corrected map indicates that this spot remains prominent, 
albeit at a level completely consistent with Gaussian statistics. Consequently, it is 
debatable how much emphasis should be placed on this residual feature, but we con- 
sider the effect of modestly increasing the scaling of the template. A factor of only 1.2 
renders the spot indistinguishable from the background level, with no noticeable im- 
pact on the results published in Jaffe et al. for the low-1 anomalies, large-scale power 
asymmetry or wavelet kurtosis. A trivial interpretation would be that the Bianchi 
template may require a small enhancement of power on scales corresponding to the 
wavelet scale of 5 degrees. 

al. 2005, Bielewicz et al. 2005, Tomita 2005, Freeman et al. 
2005), (3) Deviations from Gaussianity, in particular, detec- 
tions in wavelet space pointing to a cold spot as the source of 
non- Gaussianity (Vielva et al. 2004, Cruz et al. 2005, Cayon 
et al. 2005). 



Several works have explored the possibility of account- 
ing for some or all of the anomalies observed in the WMAP 
data by introducing some shear and vorticity in the universe 
(Jaffe et al. 2005a,b, Land & Magueijo 2005b, McEwen et 
al. 2005b). Models introducing these anisotropic character- 
istics fall under the class of Bianchi Vllh models (Barrow 
et al. 1985). In this paper we confirm the validity of the 
Bianchi model obtained as a best fit to the WMAP data 
by Jaffe et al. 2005b, as a possible explanation of all the ob- 
served anomalies. Following the same procedure as in Cayon 
et al. 2005 we show that the Higher Criticism HC statistic 
(Donoho & Jin 2004) of the Bianchi corrected WMAP map 
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is compatible with Gaussianity (at the 99% c.l.) at all the 
wavelet scales. 

This paper is organized as follows. We present the for- 
malism in Section 1. Section 2 is dedicated to introduce the 
analyzed data and simulations. Results are presented in Sec- 
tion 3 and they are discussed in Section 4 (conclusions also 
included in this section) . 



2 FORMALISM 

As indicated above we here follow the formalism presented in 
Cayon et al. 2005. The statistical study is based on the HC 
statistic proposed in Donoho & Jin 2004, Jin 2004. For a set 
of n individual observations Xi from a certain distribution, 
HC is defined as follows. The Xi observed values are first 
converted into p-values: pt = P{\N(0, 1)| > |A;|}. After 
sorting the p-values in ascending order p( X ) < p( 2 ) < • • • < 
P(„j , we define the HC at each pixel i by: 



VP(i)( 1_ P(0) 

Unusually large values of HC n ,i imply deviations from Gaus- 
sianity. The fact that the statistic is calculated at every pixel 
will allow for the location of the source of any detected devi- 
ation. In order to quantify the statistical level of any detec- 
tion the maximum of the HC n ^ for a given map (generally 
denoted HC) will be compared to the distribution of HC 
values obtained from Gaussian simulations (see discussion 
in the next section). HC has been shown to capture the 
unusual behavior of the few most extreme observations as 
well as any unusually large number of moderately high value 
observations. One should note that this is a completely dif- 
ferent test to one based on the kurtosis statistic. The kurto- 
sis enhances any deviation from Gaussianity hidden in the 
fourth order moment and it depends on the bulk of the data. 

The HC of the Bianchi corrected WMAP data is es- 
timated in real and in wavelet space. Several groups have 
presented analyses of the WMAP data in wavelet space 
(Vielva et al. 2004, Cruz ct al. 2005, Mukherjee & Wang 
2004, McEwen et al. 2005a,Liu & Zhang 2005, Cayon et 
al. 2005). In particular, the Spherical Mexican Hat (SMH) 
wavelet has been shown to be sensitive to the presence of 
non-Gaussianity in the WMAP data. Analyses based on kur- 
tosis and on HC show detections above the 99% c.l. at a 
scale of five degrees. A cold spot in the southern hemisphere 
seems to be accounting for these detections. Our aim in this 
paper is to study what is the effect on the HC results caused 
by the subtraction of a Bianchi based template from the 
WMAP data. HC is therefore calculated not only for the 
Bianchi corrected WMAP map but also for the convolution 
of this map with the SMH wavelet at fifteen scales ranging 
from 13.74 to 1050.0 arcmins. 



3 BIANCHI CORRECTED WMAP DATA AND 
SIMULATIONS 

The first year of WMAP observations is available at the 
Legacy Archive for Microwave Background Data Analysis 



(LAMBDA) website^ . For the purpose of cosmological stud- 
ies the Foreground Cleaned Intensity Maps at bands Q,V 
and W are used. A combined map (the so called WMAP 
data/map in this work) is built following Komatsu et al. 
2003. The weighted combination is given by, 

The temperature at pixel i, T(i) results from the ratio of the 
weighted sum of temperatures at pixel i at each radiometer 
divided by the sum of the weights of each radiometer at pixel 
i. The radiometers Ql, Q2, VI, V2, Wl, W2, W3 and W4 
are sequentially numbered from 3 to 10. The weights at each 
pixel, for each radiometer w r (i), are the ratio of the number 
of observations N r (i) divided by the square of the receiver 
noise dispersion cr ,r- The resulting map is downgraded from 
resolution nside — 512 to resolution nside = 256 (the total 
number of pixels being 12 x nside 2 ). 

The Bianchi template that is used in this work is the 
best fit found by Jaffe et al. 2005b. They considered Bianchi 
type Vllh models characterized by the values of Q and 
x = (h/(l — Qo)) 1 ^ 2 , where h is the scale on which the ba- 
sis vectors change orientation. The best fit was obtained for 
Q = 0.5 and x = 0.62. We present here an analysis of the 
Bianchi corrected WMAP data set obtained after subtrac- 
tion of this template (and multiples of it) from the original 
WMAP data. The final analysis is performed on a masked 
version of this map that includes the so called KpO mask (re- 
leased by the WMAP team and available at the LAMBDA 
website) . An extension of this mask is applied to the wavelet 
convolved maps as discussed in Vielva et al. 2004. 

The confidence levels of the HC statistic corresponding 
to the Gaussian assumption are drawn from 5000 simula- 
tions were the power spectrum is the one that best fits the 
WMAP, CBI and ACBAR CMB data, plus the 2dF and 
Lyman-alpha data. The simulations take also into account 
the beam transfer functions, the number of observations and 
the noise dispersion for each receiver. All these are provided 
by the WMAP team through the LAMBDA website. 



4 RESULTS 

The values of the HC statistic for the WMAP and the 
Bianchi corrected WMAP data are presented in Figure 1. 
Solid, dotted-dashed and dashed lines show the 99%, 95% 
and 68% confidence levels respectively. These are obtained 
from the 5000 Gaussian simulations described above. Since 
the HC statistic of a map is a maximum value (HC = 
maXi(HC n ,i)), the confidence regions are one-sided. As one 
can see, there is no detection of deviations from Gaussian- 
ity in the Bianchi corrected WMAP data (factor of 1.0 in 
the figure). This corroborates the results obtained by Jaffe 
ct al. 2005b. Subtraction of the Bianchi template corrects 
most of the anomalies that have been observed in the first 
year WMAP data. 

The maps of HC values for the WMAP and the Bianchi 
corrected WMAP data at a wavelet scale of 5 degrees are 
presented in Figures 2 and 3. As one can immediately see 

t http://lambda.gsfc.nasa.gov/ 
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Figure 1. Values of the HC statistic for the WMAP (stars) and the Bianchi corrected WMAP (crosses) data. The bands outlined by 
dashed, dotted-dashed and solid lines correspond to the 68%, 95% and 99% confidence regions respectively. The Bianchi template is 
multiplied by different factors before subtraction. Factors 0.3,0.5 and 1.0 are considered in the results presented in the left panel (blue, 
green and magenta crosses respectively). Results for factors 1.1,1.2 and 2.0 are presented in the right panel (magenta, blue and green 
crosses respectively). 



HC WMAP sc300p0 




171 



Figure 2. HC values of the WMAP data at the wavelet scale of five degrees. 



the detected non-Gaussianity in the WMAP data is domi- 
nated by the values of a ring of pixels at the spot centered at 
I ~ 209°, b ~ —57° (note that the ring structure is likely to 
be caused by the convolution with the wavelet). Subtraction 
of the Bianchi template drops the values of the HC statistic 
at that spot below the 68% c.l. (see Figure 1). However, the 
spot still appears as one of the most prominent regions. Even 
if the relevance of this feature is debatable we decided to 
see whether different normalizations of the subtracted tem- 
plate can influence its amplitude. A factor of 1.2 reduces 
the amplitude of the spot making it comparable with the 
surrounding values as can be seen in Figure 4. We have also 
studied whether this normalization factor affects the other 



anomalies observed in the WMAP first year data. The effect 
on the alignment between the quadrupole and the octopole 
and the planarity as denned in Oliveira-Costa et al. 2004 
is presented in Table 1. These large scale anomalies are not 
very much affected by a change in the normalization factor of 
the Bianchi template. The ratio of power between two hemi- 
spheres for a certain range of I was calculated as in Eriksen 
et al. 2004a. The probability of having a maximum power 
asymmetry ratio larger than that of a given map (either the 
WMAP or the Bianchi corrected WMAP maps) is shown in 
Table 2. This probability is slightly affected by the normal- 
ization factor. There is no much of a difference between the 
results for a factor of 1 and a factor of 1.2. However, above 
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this, an increase in the normalization factor starts worsening 
the results. This agrees with the trend observed for the HC 
statistic as can be seen in Figure 1. 
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5 DISCUSSION AND CONCLUSIONS 

Several explanations have been discussed in the literature to 
account for the different anomalies observed in the first year 
of WMAP data ( Schwartz et al. 2004, Vale 2005, Freeman 
et al. 2005, Land & Magueijo 2005b, Jaffe et al. 2005a,b, 
Tomita 2005). These anomalies include: (1) The low value 
of the quadrupole, its alignment with the octopole and other 
multipole alignments, (2) The asymmetries observed in the 
distribution of large scale power, (3) The non-Gaussianity 
detected in wavelet space which source seems to be centered 
at a cold spot in the southern hemisphere. However, much 
of this work has concentrated on explaining the anomalous 
low-1 structure. Jaffe et al, however, serendipitously discov- 
ered that at least part of the large angular scale structure of 
the WMAP data was described by a particular Bianchi Vllh 
model, and that correcting the data for this contribution 
provided a potential resolution for each of the above anoma- 
lies. We here confirm their results by estimating the effect of 
subtracting the proposed Bianchi template from the WMAP 
data, on the HC statistic. A detection of non-Gaussianity lo- 
calized in a cold spot at I ~ 209° , b ~ —57° was found by the 
estimation of the HC statistic of the WMAP data (Cayon 
et al. 2005). Subtraction of the Bianchi template renders 
the WMAP data statistically compatible with the expected 
levels based on Gaussian simulations. 

The map of HC values for the Bianchi corrected 
WMAP data (see Figure 3) still shows the cold spot in 
the southern hemisphere as one of the most prominent re- 
gions. Although no special consideration should now be 
given to the spot on a statistical basis, it remains interesting 
given that the original detection of non-Gaussianity was inti- 
mately associated with the region. We have determined that 
a modest increase in the amplitude of the Bianchi template 
(by a factor 1.2) renders the spot indistinguishable from the 
background level. Importantly, such a change in amplitude 
does not perturb the effect of the template in resolving is- 
sues related to the low-1 multipole features and power asym- 
metry. There is no compelling statistical reason, therefore, 
to believe that the Bianchi template amplitude derived in 
Jaffe et al should be considered underestimated. The HC re- 
sults presented here, however, may be interpreted as imply- 
ing that the Bianchi template may require a small enhance- 
ment of power on scales corresponding to the wavelet scale 
of 5 degrees. One could create a Bianchi template with more 
power on small scales by reducing the matter density, but 
this would also affect the large-scale structure and therefore 
the significance of the fit to the data. The most pragmatic 
interpretation of Jaffe et al - that the best-fit Bianchi model 
provides a template temperature pattern which alternative 
models would need to reproduce in order to resolve the ob- 
served anomalous anisotropy structure - remains valid. 
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Table 1. Alignment and planarity of low order multipoles. The normalization F of the Bianchi template is included in the first 
column. Second and third columns correspond to the direction of I = 2 (Galactic longitude and latitude in degrees) . The direction 
of I = 3 is shown in the fourth and fifth columns (I and b in degrees). The angle a (degrees) corresponding to the alignment between 
the quadrupole and the octopole is shown in column 6. The probability of finding a weaker alignment is presented in column 7. 
Values in the last three columns correspond to the probability of finding a more planar multipole for / = 3, 5, 6 respectively. 
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h 


f>2 


; 3 


63 


a 


P 


P 3 


Ps 


Pe 


0.0 


61.9 


247.8 


63.4 


232.8 


7.0 


0.992 


0.89 


0.001 


0.98 


0.5 


8.3 


57.1 


73.7 


253.3 


82.7 


0.127 


0.73 


0.003 


0.99 


1.0 


3.3 


44.5 


56.8 


323.5 


82.3 


0.134 


0.55 


0.05 


0.98 


1.2 


0.3 


40.6 


45.6 


327.2 


78.3 


0.203 


0.62 


0.10 


0.97 


2.0 


9.0 


212.2 


26.9 


317.3 


80.9 


0.159 


0.86 


0.37 


0.92 



Table 2. Probability of the observed asymmetry in the power spectrum between the two hemispheres. The normalization factor F 
is given in column 1. The probability for I = 2,20 is presented in column 2. Values in column 3 correspond to the probability for 
I = 2,40. 



F P(/ = 2,20)(%) P(/ = 2,40)(%) 



0.0 1.0 0.7 

0.5 9.8 6.0 

1.0 17.0 22.0 

1.2 15.9 28.9 

2.0 4.3 13.6 
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Figure 3. HC values of the Bianchi corrected WMAP data at the wavelet scale of five degrees. 

HC WMAP-1.2*Bianchi sc30CpO 




0.0 ^^^^^^^m ^^^^^h 19.5 



Figure 4. HC values of the WMAP data minus 1.2 times the Bianchi template, at the wavelet scale of five degrees. 
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